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Abstract 



The zero temperature vortex phase diagram for uniaxial anisotropic su- 
perconductors placed in an external magnetic field tilted with respect to the 
axis of anisotropy is studied for parameters typical of BSCCO and YBCO. 
The exact Gibbs free energy in the London approximation, using a self-energy 
expression with an anisotropic core cutoff, is minimized numerically, assum- 
ing only that the equilibrium vortex state is a vortex-line-lattice with a single 
vortex line per primitive unit cell. The numerical method is based on sim- 
ulated annealing and uses a fast convergent series to calculate the energy of 
interaction between vortex lines. A phase diagram with three distinct phases 
is reported and the phases are characterized in detail. New results for values 
of the applied field close to the lower critical field are reported. 
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I. INTRODUCTION 



One important theoretical problem in the study of type-II superconductors is the zero 
temperature equilibrium vortex phase diagram for clean, uniaxially anisotropic, bulk mate- 
rials subjected to an external magnetic field, H tilted with respect to the axis of anisotropy 
(c-axis). In the case of isotropic materials the vortex phase is the well known triangular 
vortex-line lattice with the vortex lines oriented parallel to H. In anisotropic materials the 
vortex phase is still a vortex-line lattice but, in general, the lattice is not triangular nor are 
the vortex lines oriented along H 

Previous work find that the vortex phases for the anisotropic superconductor are related 
to those for the isotropic one when H is oriented along the c-axis, in which case they are 
identical, when H is along the a-b plane or for H large compared with the lower critical 
field 0. In the latter two cases, the vortex-line-lattices are related to the triangular lattice 
by scaling and the vortex lines are essentially parallel to H. For general orientations of H 
this is not the case. The vortex phases were first obtained by Buzdin and Simonov 0] by 
numerically minimizing the Gibbs free-energy in the London limit with material parameters 
appropriate for YBCO. Our calculation, like their's, minimizes the London Gibbs free-energy 
numerically. However, we go beyond the work of these authors by incorporating into our 
calculations recent results for the vortex state at the lower critical field, not taken into 
account by them [^,0. We show that this leads to new vortex phases close to the lower 
critical field. 

It was found first by Grishin, Martynovich and Yampol'skii and later by others [p|,|9|,p!0 



that a pair of straight vortex lines parallel to each other and coplanar with respect to the 
c-axis attract at large distances. A consequence of this is that a chain of such lines has 
lower energy than the same vortex lines placed far apart. On the basis of this result Buzdin 
and co-workers P,p!0| suggested that at the lower critical field the mixed state consists of a 
vanishing small density of vortex line chains instead of vortex lines. However this argument 
does not apply in general. As pointed out in Ref. 0, in order to determine whether vortex 
lines or vortex-line chains constitute the vortex phase at the lower critical field it is necessary 
to minimize the Gibbs free-energy for each one of them and obtain their respective lower 
critical fields. The one with the lower critical field is the true equilibrium phase. One 
important aspect of the calculation of Ref. |^ is the use of an expression for the tilted 
vortex-line self energy with an appropriate anisotropic core cutoff function. This expression 
was derived by Subdo and Brandt |11| and later used in Ref. ^ to show, by minimizing the 
London Gibbs free-energy for a vanishing small density of vortex lines, that for anisotropy 
parameters typical of BSCCO there is coexistence two such states, with different vortex lines 
tilt angles, for one particular orientation of H. These authors also showed that if a self- 
energy expression based on an isotropic core cutoff, such as the one adopted by Buzdin and 
Simonov is used instead, the coexistence disappears. It is found in Ref. |^ that at the 
lower critical field and for anisotropy parameters typical of BSCCO the equilibrium vortex 
state for H tilted hj a < ac = 7.9° with respect to the c-axis is a dilute vortex-line chain 
state whereas for a > etc = 7.9° it is a dilute vortex-line one. For anisotropy parameters 
typical of YBCO Ref. |^ finds that the lower critical fields for vortex lines and vortex-line 
chains differ so little from each other that it is not possible do decide which one corresponds 
to true equilibrium. 
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In this paper, we minimize the exact Gibbs free-energy in the London hmit using Sudbo 
and Brandt [11| expression the vortex hne self energy. We assume that the equihbrium state 
the vortex hues are straight and parallel to each other and the VLL has only one vortex line 
per unit cell, and obtain the vortex lines orientation and VLL shape as functions of H, by 
minimizing the Gibbs free-energy in the London limit. 

The free parameters in our calculation for given H, specified here by its magnitude, H, 
and its tilt angle with respect to the c-axis, a, are the vortex lines tilt angle, 9, and the 
VLL primitive unit cell vectors Li and L2, defined in the plane perpendicular to the vortex 
lines. Since the vortex lines are straight and parallel to each other, the magnetic induction 
B is also tilted by 9 with respect to the c-axis and, as explained in detail in Sec. |11 A|, 



it is coplanar with H and the c-axis. The orientations of H, B, and the c-axis, and the 
VLL primitive unit cell are shown in Fig. |1|. Hereafter we specify the VLL unit cell by the 
magnitudes of Li and L2, denoted by Li and L2, by the angle between Li and L2, v^, and 
by and the angle of rotation of the VLL unit cell with respect to B, ip, as shown in Fig. ^ 
As discussed in Sec. [II A| , a minimum of the Gibbs free-energy occurs for = 0, where if) is 
defined in Fig. |1|. We assume that this is the absolute minimum and carry out the calculations 
with ?/' = 0. This means that Li is also coplanar with H and the c-axis,oriented along the 
x-direction defined in Fig. |I]. To simplify the calculations we set 99 = 90°. This means that 
the VLL unit cell is a rectangle, with L2 perpendicular to the plane defined by H, B and the 
c-axis (Fig. |l|). We find that this simplifications does not change our main results. Under 
these conditions we have only three free parameters, to be determined by minimizing the 
Gibbs free-energy: 9, Li and L2. This minimization is carried out numerically, using a fast 
convergent series to calculate lattice sums entering the expression for the VLL interaction 



energy [|T3|. This series, together with the small number of free parameters, allow us to 
create an efficient and accurate numerical algorithm. 

We study in detail the phase diagram for anisotropy parameters typical of YBCO and 
BSCCO. For YBCCO our results are essentially the same as those obtained by Buzdin and 
Simonov, except very close to the lower critical field. For BSCCO our results are new. 

The main results of this paper are summarized in the generic phase diagram show in 
Fig. 1^. We find three distinct phases, hereafter called phase-I, phase-II, and phase-Ill. We 
find that phase-I is new and closely related to the vortex-line equilibrium state at the lower 
critical field. We also find that phase- II is essentially identical to that found by Buzdin 
and Simonov and that phase-Ill is that expected for fields much larger than the lower 
critical field. In this limit the VLL unit cell is related to that appropriate for the isotropic 
superconductor by scaling, and B ~ H (^^ ~ a) |^. 

The lower critical field in isotropic superconductors is independent of the orientation 
(a), and has magnitude Hd- For H just above ifci, the vortex state corresponds to vortex 
lines parallel to H very far apart. In our notation this state is described by 6' = a and Li, 
L2 = 00. In a superconductor with uniaxial anisotropy, the magnitude of the lower critical 
field depends on a, hereafter denoted by Hci{a) f^. The vortex state for if just above Hci{a) 
consists of vortex lines not parallel to H (^^ 7^ a) in one of two possible configurations, named 
dilute vortex line state (DVL) and dilute vortex-line chain state (DVLC) 0. The DVL is 
similar to the traditional vortex state for the isotropic superconductor, in the sense that Li, 
L2 = 00, whereas the DVLC corresponds to vortex-line chains along the x-direction with 
period Li, very far apart, that is Li finite and L2 = 00. The latter possibility results from 
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attractive interactions between vortex lines that only exist in anisotropic superconductors 
i7|j8|J^,p!0f| . To decide if the DVL or the DVLC is the equilibrium vortex state it is necessary 
to calculate Hci{a) for each one. The equilibrium configuration is that with the smallest 
Hci{a). This calculation is reported in Ref. 0. The results for anisotropy parameters 
typical of BSCCO and YBCO, referred to as strong an moderate anisotropies, respectively, 
are shown in Fig. |^. It is found that for both anisotropies there is a value of the external field 
tilt angle, ac, around which the equilibrium vortex states undergoes a significant change. 
For a > Oc and both for strong and moderate anisotropies the equilibrium state is a DVL 
with vortex lines tilted at ^ ~ 90°, that is parallel to the ab-plane. For strong anisotropy 
the vortex state for a < Oc = 7.9° is a DVLC with vortex lines tilted at 9 < 60°. At 
a = ac, there is coexistence of a DVL and DVLC. A discontinuous change in 6 takes place 
at a = ac- For a moderate anisotropy, 6 is found to change smoothly with a and ac ~ 20°. 
The calculation on Ref. |0 finds that for moderate anisotropy the DVL and DVLC values 
of Hci{a) are so close that their difference is within the numerical errors of the calculation. 

We find that phases I and II evolve, respectively, from the a < ac and a > ac vortex 
states at H = Hci{a), represented by the horizontal line in Fig. ^. 

Phase-II is essentially an extension into the H > Hci{a) region of the phase diagram of 
the vortex state found for a > ac at H = Hci{a). Its predominant characteristic is that 
6 ~ 90° and, consequently, the VLL within phase-II is essentially identical to that for vortex 
lines parallel to the ab-plane. In this case the unit cell is, approximately, related to that for 
the isotropic superconductor by scaling relations that depend only on the anisotropy. For 
strong and moderate anisotropies the unit cell is a rectangle elongated in the y-direction, 
that is Li ^ L2. Phase-II is limited from above by the line H = H*{a), which gives the 
field value around which 6 changes from 90° to a lower value, accompanied by significant 
changes in the VLL unit cell. At a = ac, H*{a) = Hci{a). 

Phase-I is the region of the phase diagram where the VLL properties are strongly in- 
fluenced by the attractive interactions between pairs of vortex lines aligned along the x- 
directions. For both strong and moderate anisotropies we find that phase-I is not limited to 
a < ac, but it also exists in a small part of the region a > ac and H > H*{a). For strong 
anisotropy we find that phase-I consists of vortex-line chains along the x-direction, with a 
period Li that is nearly identical to that of an isolated chain with the same tilt angle (which 
differs from that of an isolated vortex), separated from each other by L2 ^ Li. This means 
that the attractive interactions between vortex lines prevail in determining the VLL period- 
icity along the x-direction. We also find that as H ^ Hci{a) from above, the equilibrium 
VLL approaches the same DVLC found in Ref. for = Hci{a). For moderate anisotropy 
we obtain a new result for the vortex state at H = Hci{a). By extrapolating our data as 
H —>■ H*{a) from above, we conclude that the vortex state at H = Hci{a) is a DVLC. 
However, within phase-I for moderate anisotropy, the attractive interactions between vortex 
lines are so weak that vortex-line chains along the x-direction with periods close to those of 
isolated chains only exist within phase-I near H = Hci{a). However, distinctive behavior, 
due to attractive vortex-line interactions, can still be identified within phase-I. 

This paper is organized as follows. In Sec. [II A| we review London theory for the energy 
of the vortex-line system under consideration and the known results for the phase diagram 
at H = Hci{a) and for particular cases at H > Hci{a). In Sec. [II B| we discuss the nu- 
merical method used to minimize the Gibbs free energy. In Sec. ^ we report the results 
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of this minimization, the detailed properties of each phase, and discuss how the phases are 
identified. Our conclusions are stated in Sec. IV. In the Appendix the relevant formulas for 
the fast convergent series used in our in our calculation are summarized. 



II. METHOD OF CALCULATION 

Here, we review London theory for this problem, its known results for the phase diagram, 
and describe some details of our numerical method. 



A. London Theory 

To obtain the zero-temperature phase diagram for a superconductor in an applied field 
it is necessary to determine the vortex arrangement that minimizes the Gibbs free-energy 
density 

with H and the volume kept constant. In Eq. ([1|) F is the vortex energy density and B is the 
magnetic induction. We assume that the zero-temperature equilibrium vortex configuration 
consist of straight vortex lines, parallel to each other, arranged on a periodic VLL, with a 
single vortex line per primitive unit cell. In this case B is parallel to the vortex lines direction 
and has magnitude B = ^q/Ac, where Ac is the VLL primitive unit cell area. Minimization of 
G gives then the vortex lines orientation and unit cell shape that corresponds to equilibrium. 
Because F is unchanged, if the vortex lines are rotated around the c-axis, the equilibrium 
B lies in the same plane as H and the c-axis. Thus, for H with magnitude H and tilted 
with respect to the c-axis by an angle a, G depends on the vortex lines tilt angle with 
respect to the c-axis, 6, and on the parameters defining the VLL primitive unit cell. In the 
system of axis defined in Fig. [^.a, the generic VLL primitive unit cell lies in the x — y plane, 
and is defined by the primitive unit vectors Li and L2, as shown in Fig. |^.b. In this case 
Ac = L1L2 sin (f. 

In the London limit the superconductor is characterized by the penetration depths Xab 
and Ac, for currents parallel to the a6-plane and to the c-direction, respectively. The free 
energy density of a generic arrangement of these vortex lines can be written as the sum of 
pairwise interactions |0] 

(2) 

where Tj is the i-th vortex-line position vector in x — y plane, A is the sample area in this 
plane, and /(r) is the Fourier transform of 

1 -^ \^k^ 

' (l + A^,P)(l + AiA:2 + A2A;2) ' 

where 
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= A^;, sin^ 9 + Xl cos^ 9 , (4) 
and gik) is the vortex core cutoff function It is costumery to write G as 

G{e, Li, L2; H) = ^[£sf (^) + eUO, Li, L2) - cos {6 - a)] , (5) 

where £sf(^) is the vortex hne self-energy and ein(6', Li, L2) is the interaction energy per 
vortex hne. In terms of /(r), defined in Eq. these quantities are given by 

esAe) = ^f{r = 0), (6) 



and 



$0 



where Rj denotes the VLL positions (Rj = niLi + n2L2, ui, n2 = integer). 
The minima of G, Eq. (J^), are known in the following cases. 

i) Lower critical field: H = Hci{a). The equilibrium vortex phase at the lower critical 
field is studied by Sudb0, Brandt and Huse 0, assuming that it consists of a dilute ar- 
rangement of straight vortex lines (DVL). Accordingly, these authors minimize G, Eq. 
neglecting ^in, and using for ^sf the expression derived by Sudb0, and Brandt [O], based on 
an elliptic core cutoff function, 

/ / N , A^COS^e 1^ i>^l + >^'e) ^ ro\ 

Ue) M^h) + ^^j^^jj^ In ] , (8) 

where Eq = (<l>o/47rAafe)^, 7 = ^ab/K = idiab, 1^ = Kb Hah, iab being the a6-plane coherence 
length. This line energy is obtained from the anisotropic Ginzburg-Landau theory using 
the Klemm-Clem transformation [jlT],|l2|. The elliptical core cutoff has semi major axis 
and semi minor axis with = cos^ ^ + sin^ 0. Sudb0, Brandt and Huse obtain 
Hci{a) for several values of the anisotropy parameter 7 and of k, and find that for certain 
ranges of 7 and k, there is a value of a, a = Oc, for which coexistence of two DVL differing 
from each other by the equilibrium 9 takes place. 

In a recent publication by some of us ||^, the Sudb0, Brandt and Huse calculation 
was generalized to account for two possible vortex phases at the lower critical field: a dilute 
vortex-line state (DVL) considered by them and a dilute vortex-line chains arrangement 
(DVLC). The DVLC consists of vortex lines tilted with respect to the c-axis, parallel to one 
another, and aligned along the x-direction, forming a periodic chain. This is possible because 
the interaction between a pair of such vortices is attractive at large distances In 
Ref: G, Eq. (|]), is minimized in the limit of vanishing vortex density, for each one of 
these possibilities. For the DVL, ^int is neglected in Eq. (|), whereas for the DVLC ^int in Eq. 
(§) is identified with the interaction energy per vortex line of an isolated chain. The lower 
critical field as a function of the tilt angle a, H(.i{a) is then obtained for each possibility, and 
the equilibrium vortex phase aX H = Hci{a) is identified as that with the smallest H(.i{a). 

In Ref: 0, two sets of anisotropy parameters are considered: 
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Strong anisotropy: k = 10, 7 = 1/ -\/200 
Moderate anisotropy: k, = 50, 7 = 1/5 

These values for 7 are typical of YBCO (moderate anisotropy) and BSCCO (strong 
anisotropy). The parameter k for moderate anisotropy is typical of YBCO but for strong 
anisotropy is about five times smaller than those typical of BSCCO [|1|. However, this dif- 
ference does not alter significantly the phase diagram because only the self-energy depends 
on K. 

The main results of Ref. [^, summarized in Fig. |^ are as follows. For both moderate 
and strong anisotropy and for a > Oc the vortex lines are nearly parallel to the a-b plane 
{6 ~ 90°). For strong anisotropy Oc = 7.9° and for moderate anisotropy a ~ 20°. For strong 
anisotropy and a < ac the vortex phase aX H = Hci{a) is a DVLC. For a > ac, on the 
other hand, no significant difference is found between the DVLC and DVL phases, because 
the vortex lines are tilted at ^ ~ 90°, and at such tilt angles the equilibrium chain period 
is very large. Phase coexistence is also found for strong anisotropy at a = = 7.9°, but 
the coexisting phases are a DVLC and DVL, instead of the two DVL found in Ref. [^]. A 
discontinuous jump in 9 takes place as a crosses ac = 7.9°. For moderate anisotropy the 
DVL and DVLC phases give nearly identical values for Hci{a), so that no conclusion on 
which is the equilibrium state is reached. 

ii) a = 0° and a = 90°. In both cases 6 = a and the lower critical field is related 
to the self energy by the usual formula Hci{a = 0°,90°) = 47r£:sf(6' = 0°,90°)/$o. For 
a = ^ = 0° the equilibrium VLL is the familiar isotropic superconductor triangular lattice, 
with Li = L2 = La = (<I>o/5 sin 60°)^/^ Lf = 60° and ip undetermined. For a = 9 = 90° 
the x-direction coincides with the negative c-axis. The equilibrium VLL is related to that 
for q; = ^ = 0° (assuming Li along the x-direction or ip = 0°) by the following scaling 
relations || 

Li = Lo7^/^x 

L2 = ivo[7^^^ cos V2X + J'^^"^ sin ipy] , (9) 
where Lq = {^0/ B siiK^Y^'^ , with Lp = 60° for the triangular lattice. These results also 



apply if the a = 0° VLL is a square lattice, which will be considered in Sec. |ITI| , in which 
case = 90° in Eq. (^. The unit cell defined by Eqs. (^ results from compressing the 
a = ^ = 0° VLL unit cell (with Li along the negative c-direction) by 7^/^ along the c- 
direction and stretching it by 7"^/^ in the ab plane. This transformation conserves the unit 
cell area Ac = Lq sin (p. The B vs. H curves for ^ = a = 0° and 6 = a = 90° are also related 
by scaling as follows |]T],0 

g-ga(90°) H-Hc^m . 

iii) High fields {H ^ Hci{a)). To a good degree of approximation the equilibrium VLL 
is related to the a = 6* = 0° VLL by the scaling relations Ref. 

U = lA^i^ X 

Ac 

L2 = Lo[(^)^/^cosy. x+ (^)V2sin¥, y] . (11) 

Ac At) 
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In this case B ~ H, with B = $o/-^a sin with (p = 60° for the triangular VLL and (p = 90° 
for the square VLL. 

Apart from these cases, httle else is known about the vortex phase diagram. Results for 
H = Hci{a), such as the occurrence of a DVLC as the equilibrium state, phase coexistence, 
and large differences between 9 and a, suggest a low field {H > Hci{a)) phase diagram that 
is non-trivial and that differs considerably from the high-field one. It is of interest then to 
obtain a full picture of the phase diagram. This is our main motivation here. 

According to the discussion above, for a given H, G depends on five independent vari- 
ables, 6',Li,L2, or, equivalently, on 9, Li, L2,(p,ip (Fig- Q-^)- It can be shown that the 
dependence of G on the angle of rotation of the VLL with respect to B, ip, is such that G 
has a minimum for ip = 0. We assume that this is the absolute minimum and, from here on, 
take ip = in our calculations. In order to simplify the calculations, we fix the value of the 
angle (f at ip = 90°. Thus, in the results reported here, Li is along the x-axis and coplanar 
with B, H and the c-axis, L2 is along the y-axis, and the VLL unit cell is rectangular. We 
find that setting ip = 90° does not change the general picture for the phase diagram. With 
these restrictions G depends only on three variables: 9, Li, L2. In the remainder of the paper 
we use instead 9,Li,B = $o/(-^i-^2) as independent variables. 

One novel aspect of our calculation is to account for the effects of interactions between 
vortex lines in the zero-temperature phase diagram. We find that the simple model described 
above has a non-trivial low-field phase diagram because attractive interactions between 
vortex lines give an important contribution to ^int at low fields. The reason is that, as 
discussed above, the equilibrium VLL unit cell has Li along the x-direction {ip = 0). This 
means that the equilibrium phase consists of periodic vortex-line chains aligned along x, 
with period Li, separated from each other by L2. According to []7|,|8|,|9|JT0|] , the interaction 
between a pair of vortex lines in the same chain is attractive if their separation x is greater 
than Xm{9), and repulsive if x < Xm{9). Thus, as long as Li is not too small compared to 
Xm{9), the vortex-line chains contribution to ^int is predominantly attractive and essentially 
negative. Interactions between interchain vortex lines are always repulsive, thus giving a 
positive contribution to ^int- 

In the DVLC state at if = Hci{a), these chains are infinitely far apart (L2 = 00) and 
Li minimizes ^int for the equilibrium 9. In this paper, we obtain the vortex phases at 
H = Hci{a) by extrapolating our results for H > Hci{a). For strong anisotropy, our results 
are in agreement with those of Ref. 0. For moderate anisotropy, we find that the vortex 
phase at if = Hci{a) is also a DVLC for a < Uc, a result not obtained in Ref. 

Our results for the zero-temperature phase diagram can be summarized in the generic 
phase diagram show in Fig. |^, consisting of three distinct phases, whose main characteristics 
are discussed in Sec. |. 



B. Numerics 

The main difficulty to carry out the minimization of G, Eq. (^, is to evaluate the lattice 
sum in the expression for e-m, Eq. (|^). This can only be done numerically. However, the 
direct evaluation of this sum is impractical for our purposes, because it converges very slowly. 



To circumvent this difficulty we use the results obtained by Doria |13|, who showed that the 
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lattice sum in Eq. can be transformed into a series that converges much faster than the 
direct sum. 

Our numerical method uses a simulated annealing algorithm to locate the minima of G. 
For a given H, specified by H and a, the minimization of G is carried out in the space 
spanned by 6', Li, and B. Our algorithm starts from a convenient choice for these variables 
and attempts to change them to new values that decrease G. Once a set of such values is 
found, another one is searched by repeating the procedure, and so on. This converges to 
a set corresponding to a minimum of G after several steps. For each new set of d,Li,B 
we calculate £in, using Doria's fast convergent series. This allows us to run the simulated 
annealing for a very large number of steps, and to obtain the minima of G with high accuracy, 
except when H is very close to Hci{a). In this case our method converges slowly, because 
the VLL unit cell becomes very large and Ei^ very small. A similar numerical method is 



described in Ref. 14 



III. RESULTS 

Here we report in detail the results obtained by the numerical method described in Sec. 



11 B| for strong and moderate anisotropics. To explore the {H, a) phase diagram we fix a 
and obtain the minima of G for several values oi H > Hci{a). By repeating this procedure 
for several a we obtain the phase diagram. We check that these minima are absolute ones 
by running the above described simulated annealing algorithm starting from distinct initial 
states, and selecting the minima with the smallest G. 



A. Strong anisotropy 

The main results are summarized in Figs. ^ and |^. 

First, we argue that our results are consistent with the prediction, discussed in Sec. |, 
that at if = Hci{a) the vortex configuration for a < Oc = 7.9° is a DVLC 

The existence of a DVLC ai H = Hci{a) for a < ac means that as H/Hci{a) 1 the 
primitive unit cell sides behave as Li Lch{= chain period), L2 00, and that Em, Eq. (^, 
approaches a negative value, equal to the interaction energy per vortex line in an isolated 
chain. Our results for a = 7° show exactly this behavior. The L1/L2 curve in Fig. ^.c 
decreases continuously as H/Hci{a) — >• 1, indicating that L^jLi extrapolates to L1/L2 = 
ai H = Hci{a). On the other hand, the Li curve in (Fig. ^.b) approaches a finite value . We 
find that the extrapolated values of Li, 9 and Em at H = Hci{a) are Li/Xab = 2.8, 9 = 50° 
and Ein/Eo = —0.04. These values agree, within numerical errors, with the chain period, 
tilt angle and interaction energy per vortex line of the isolated vortex-line chain obtained in 
Ref. §. 

Next we consider the behavior for a < ac and H > Hci{a). 

For a = 7° our results for the Li vs. H curve (Fig. |^.b) show that there are two distinct 
behaviors. One for 1 < H/Hci{7°) < 1.03, where Li increases with H, and another for 
H/Hci{7°) > 1.03 where Li decreases with H. The vs. H curve changes slope at 
H/Hci{7°) = 1.03, and shows that Em increases continuously with H, being negative for 
H/Hci{7°) < 1.02. A decrease of Li with H is what is expected if interactions between 
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vortex lines are purely repulsive. This behavior is identified with phase-Ill (Fig. |]) and 
corresponds to the region where L1/L2 > 0.6 , 9 < 23° and > 0. The increase of Li 
with H, on the other hand, arises from attractive interactions between vortex lines. This 
can be seen by plotting Li vs. 6 for a = 7°, as shown in Fig. ^.c. This figure shows that 
in this range of H values, Li is nearly equal to the period of an isolated vortex-line chain 
with the same 6 and has essentially the same H dependence. This is only possible if the 
attractive intrachain vortex-line interactions are strong compared to interchain repulsions. 
Further evidence for this is found in the negative values for Sin- We identify the region 
1 < H/H^i{7°) < 1.03 as phase-I and that for H/Hci{7") > 1.03 as phase-Ill. For a = 7° we 
also find that L2 > Li, and that L2 decreases with H in such a way that L1/L2 increases 
with H as shown in Fig. ^c. 

A similar behavior is expected for other values of a < ac- For a = 5" we find that the 
L1/L2 curve shows a downturn as H/Hci{a = 5°) decreases (Fig. §.c), similar to that for 
a = 7°, suggesting that the limit value as H/Hci{a) 1 is small. We also find that the limit 
value of 9 as H/Hci{a) 1 (Fig. |^.a) agrees well with that obtained in Ref. 0. Although, 
for a < 5°, we are unable to observe a region where Li increases with H and Em < 0, we 
believe that the phase diagram for a = 5° and a = 7° are similar. One difficulty with values 
of a close to 0° is that the region where phase-I exists is so close to H/Hci{a) = 1 that the 
accuracy of our numerical method is insufficient to observe it. 

Next we discuss the behavior for a > ac and H > Hci{a). 

Our results show that 6 remains close to 90° in the range Hci{<y) < H < H*{a) (Fig. 
§.a). We find that for the results in this range the scaling relations discussed in Sec. |[ For 
strong anisotropy 7 = l/v200, so that tan^ > I/7 for 6 > 86°. Accordingly, the field 
H*{a) is obtained from the 6 vs. H curve as the largest H for which 6 > 86°, as indicated 
in Fig. g.a. For a > 10° and Hci{a) < H < H*{a) the L1/L2 and B vs. H curves in Fig. 
§ essentially coincide. The common curves are L1/L2 = 7 = 0.07 and the B vs. H curve 
for a = 9 = 90°, obtained from that for a = = 0° by the scaling relations Eq. (jT^). The 
behavior of Li, L2 and 6^ in this region is that expected for repulsive interactions: as H 
increases, Li and L2 decrease, and Bin increase (Fig. ^). From these results we identify the 
region -ffci(a) < H < H*{a) as phase-II. 

For H > H*{a) and for a = 10°, 15° and 20° the tilt angle 6 decreases smoothly with 
increasing H (Fig. ^.a). The behaviors of Li and Em with H are not smooth. As shown in 
Fig. ^ there are regions where Li increases with H and the Ein curve shows a dip, reaching 
negative values for a = 10°. This behavior of Li is similar to that for a = 7° < in 
phase-I region. We also find that for a < 20° the VLL period along x, Li, is close that for 
an isolated vortex-line chain with the same 6, and has a similar H dependence (Fig. |^.c). In 
view of these similarities we identify the region where Li increase with H as phase-I. Thus 
the portion of the phase diagram occupied by phase-I extends into the region a > ac for 
H > H*{a), as shown in Fig. ^ 

Our results indicate that for strong anisotropy within phase-I the vortex-line chains along 
X have essentially the same period as an isolated chain tilted at the equilibrium 6. 

We find that as a increases above a = 20° the region where Li increases with H becomes 
smaller and the dip in Ein decreases. Eventually these features disappear altogether at some 
a, so that the phase-I region in the phase diagram is bounded as shown in Fig. |^. 

Our results indicate that the transition from phase-II to phase-I is a smooth crossover, 
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with an intermediate region starting just above H*{a) and ending when 6 ~ 70°. In this 
region the L1/L2 and B curves for a = 10°, 15° and 20° show only small deviations from 
the scaling behavior of phase-II, with L1/L2 smaller than the scaling value and the B curve 
slightly above the scaling one (Fig. Behavior characteristic of phase-I, as described above, 
shows up only for 6 ^ 70°. The reason is is that for 6 ^ 70° interactions between intrachain 
vortex lines are predominantly repulsive because Li <^ Xm{G)- 

For a = 10°, 15° and 20° phase-I behavior ends at the H value where Li stops increasing 
with H. For larger H, we find that Li decreases smoothly with H, with L1/L2 ~ 1, and 
that Bin increases with H (Fig. ^. We identify these regions as phase-Ill (Fig. ^. 

B. Moderate anisotropy 

The results are similar to those for strong anisotropy and are summarized in Figs. ^ and 

First we argue that as H/Hci{a) 1 the equilibrium state approaches a DVLC for 
a ^ 20°. As seen in Fig. ^.c, for a = 10°, 15° and 20°, L1/L2 approaches a small value as 
H/Hci{a) — s> 1, which we interpret as indicating that L2 00. The Li vs H and Li vs. 9 
curves, shown in Figs. |^.b and 0.c, suggest that, as H/Hci{a) 1, Li approaches a finite 
value that is close to the period of an isolated vortex-line chain with the same 6 ^ . The 
vs. H curve shown in Fig. |^.a becomes negative as H/Hci{a) — >■ 1 for a = 10° and 13° and 
extrapolates to a negative value at H = Hci{a) for a = 15° and 20°. This behavior is only 
possible for a DVLC state at if = Hci{a). 

Next we consider the behavior for a < Oc and H > Hci{a). 

For a = 10°, 13° and 15° the Li vs. H curves in Fig. ^.b show non-monotonic behavior. 
As H increases Li first decreases, then increases and, finally, decreases monotonically. This 
behavior suggests that attractive interactions between intrachain vortex lines are competing 
with repulsive interchain ones. This is further supported by the negative Ein, for the same 
a values, shown in Fig. ^.b. We interpret as phase-I the H range starting from H = Hci{a) 
and ending at the H value in Fig. ^.b where Li stops increasing. However, the behavior in 
phase-I for moderate anisotropy differs from that for strong anisotropy. The Li vs. 6 curve 
in Fig. |^.c shows that Li differs considerably from the equilibrium period of an isolated 
chain with the same 6 and has a different H dependence. 

We find that for a > ac our results agree with phase-II behavior for Hci{a) < H < H*{a). 
For moderate anisotropy 7 = 1/5, and tan 6' > I/7 for 6 > 80°. The field H*{a) is obtained 
from the 6 vs. H curve as the largest H for which 6 > 80°, as indicated in Fig. ^.a. We find 
that for Hci{a) < H < H*{a) the L1/L2 and B vs. H curves shown in Fig. ^ agree well 
with the scaling predictions Sec. |. 

Thus we conclude that for moderate anisotropy the phase diagram is also like that shown 
in Fig. |. 

IV. CONCLUSIONS 

In conclusion, we obtain the zero-temperature phase diagram for superconductors with 
anisotropy parameters typical of BSCCO and YBCO. Three distict phases are found and 
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characterized in detail. Phase-I is where we find new results. The most significant ones are: 
i)in the case of strong anisotropy we determine how the vortex-line-chain state that occurs 
aX H = Hci{a) {a < etc) is modified by interchain interactions, ii) In the case of moderate 
anisotropy we show that the vortex state at if = Hci{a) is a vortex-line chain state. 

In phase-II, we recover the results of Buzdin and Simonov that reveal a large region 
of the phase diagram over which the vortex lines remain nearly parallel to the ab-plane. 
Moreover we find the curve, H = H*{a), around which the vortex lines tilt away from the 
ab-plane. We also find that this line starts at a < and H = ifci(ttc) and that, as this 
line is being crossed, the VLL undergoes a structural change that becomes more and more 
abrupt as a ac- 

Our results are obtained with the simplifying assumption that (f = 90°. We have also 
carried out the minimization of G, Eq. (|^), with (f as independent variable. The results are 
identical to those reported above, except for the VLL unit cell shape which, in general, is 
not rectangular. 

Our present results seem to indicate that the equilibrium VLL changes smoothly with H 
and a. In this case the phase boundaries in Fig. ^ are crossover lines. However, our data for 
strong anisotropy shows abrupt changes in the slope of the curves for 6, B and Li vs. H for 
H close to ifci(tt)- These could result from a rapid crossover or from true discontinuities. 
It is found in Ref. p that, for strong anisotropy, if a is changed and H is kept along the 
H/Hci{a) = 1 curve, there is a discontinuous change in 9 at a = ac = 7.9°, associated 
with coexistence of a DVLC state with a DVL one. This coexistence might remain for 
H/Hci{a) > 1 in the vicinity of ac- In the present study we do not attempt to verify if there 
are true discontinuities. However, our numerical method is capable of doing so. Work on 
this topic is under way and will be reported elsewhere. 
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APPENDIX 



In this appendix the formulas that we use to calculate ein{0, Li, L2,ip) numerically are 
summarized. 

Doria showed that the lattice sum in Eq. (0) (with Li along x) can be transformed in a 
fast convergent series [|l^. The results are: 



+ 



In 



{^Vo{0) - (1 - cos' 9 f 

Ac Jo 
]} 



du 



I cos 6*1 ,1-^61-1- Acl cos 6*1 



CiiU C2{U 



dVoju) 



+ 



X Aafe(l + I cos6'|) 
where Vq and its derivative are given by 



(Al) 



12 



^ 2/iotanh(/ioa/2) ^ 



oo . ^ 

+ 5: ( ^- -) + In (/.o/2) + (A2) 

m=l d liQ + m? ™ 



4/io ^anh (/ioa/2) sinh' (/iofx/2)^ 

2 cos (mxs) , rr, , , 1 , 

+ /^o L „2^r»2 exp(-^/i^ + m2ag)[ = + as] + 



where i?c = 0.5772... is Euler's constant and x is defined as x = {'2TrL2/ Li) cosip. The 
functions of Ci, C2, a and /iq are defined as: 

ci{u) = t| - u{l - 7^) cos^ 6 , 
C2{u) = 1-m(1-72) , 



L2 

a = 27T 



Ci{u) . 

Li \ C2[u) 



and 



2-nXc^Jciiu) 



In the expression for Eq. Vo(0) stands for Vq{u), Eq. ([A^), for -u = 0, or, in 

other words with the functions Ci, C2, a and /iq in this equation evaluated at -u = 0. 

In our numerical calculations of Bin we evaluate Vo(0) by using the following generic 
approximation for the summations in eqs.( |A^) and ([A3|) 



where g{x) is an analytic function. By using M = 200 we obtain excellent numerical precision 
with modest computer run times. 

The integral in eq . (|Al|) is calculated using standard integration algorithms found in the 
literature, up to a precision of one part in 10^. We run typical calculations of 10-20,000 
Metropolis steps for each run. The results obtained are easily reproducible. 
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FIGURES 




FIG. 1. Definition of coordinate system and VLL primitive unit cell, a) Relative orientations 
in the plane containing H, the c-axis and the equilibrium induction B. The x-axis is perpendicular 
to B and the y-axis is perpendicular to the plane H-B, pointing inside, b) Generic VLL primitive 
unit cell in the x — y plane. 
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FIG. 3. Results of Ref.[5] for Hci{a) and 6 vs. a for: a) strong anisotropy, b) moderate 
anisotropy 
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FIG. 4. Numerical results for strong anisotropy. a) Vortex lines tilt angle. The full line 9 = 86° 
corresponds to tan^ = 1/7- The arrow indicates H*{a = 10°). b) Magnetic induction, c) Ratio 
between VLL unit cell sides. The full line L1/L2 = 7 corresponds to the value predicted by scaling 
for a = 90° 
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FIG. 5. Numerical results for strong anisotropy. a) Interaction energy per vortex line. Arrow 
indicate the value of H* for a = 10". b) VLL period along x. c) VLL period along x vs. tilt angle 
compared with the period of an isolated vortex-line chain tilted at the same angle (full curve) 
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FIG. 6. Numerical results for moderate anisotropy. a) Vortex lines tilt angle. The full line 
6 = 80° corresponds to tan^ = 1/7- b) Magnetic induction, c) Ratio between VLL unit cell sides. 
The full line L1/L2 = 7 corresponds to the value predicted by scaling for a = 90° 
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FIG. 7. Numerical results for moderate anisotropy. a) Interaction energy per vortex line, b) 
VLL period along x. c) VLL period along x vs. tilt angle compared with the period of an isolated 
vortex-line chain tilted at the same angle (full curve) 
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